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Abstract. We present a 1.2 mm (250 GHz) map obtained with MAMBO on the IRAM 30m telescope of the 
central 25 arcmin^ of the proto-cluster surrounding the 2=4.1 radio galaxy TN J1338— 1942. The map reaches 
a la sensitivity of 0.6 mjy in the central area, increasing to 1.2 mjy at the edges. We detect 10 candidate mm 
sources, of which 8 are also detected in a deep VLA 1.4 GHz map and/or a VLT _R— band image. Three sources 
have a flux density 5*1. 2mm >4.0 mJy, representing a 7a overdensity compared to random field surveys, which 
predict only 1 such source in our map area. We obtained SCUBA/JCMT 850 fim and 450 fim photometry of six 
radio/optically identified MAMBO sources, confirming 5 of them with S/N>4. Radio-to-mm and mm-to-submm 
redshift estimators cannot put strong constraints on the redshifts of these MAMBO sources, but 9 of them are 
consistent within the uncertainties (mean Az—+2.6) with z—A.l. One faint MAMBO source is possibly identified 
with an extremely red object {R — K—6.1) at a likely spectroscopic redshift 2=1.18. 

The four brightest MAMBO sources are all located north of the radio galaxy, while the densest area of companion 
Lya excess and Lyman break galaxies is to the southeast. None of the 14 spectroscopically confirmed Lya emitters 
in the MAMBO field are detected at 1.2 mm; their average 1.2 mm flux density is (5'i.2mm)=0.25±0.24 mJy. If the 
mm sources lie at 2=4.1, none of them show excess Lya emission in our narrow-band images. Both populations 
thus show no apparent overlap, possibly due to dust quenching the Lya emission. If the mm sources are part 
of the proto-cluster, our results suggest that galaxies with star formation rates of a few 1000 M0yr~^ could be 
spread throughout the proto-cluster over projected scales of at least 2 Mpc. 
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* Based on observations obtained with the IRAM 30m, the 
Very Large Array, the James Clerk Maxwell Telescope, and 
the ESO Very Large Telescope at Paranal, Chile (programs 
LP167.A-0409, 69.B-0078 and 71.A-0495). 



1. Introduction 

Our knowledge of large-scale structures at high red- 
shifts has made significant progress in recent years 
with the detection of several overdense regions out to 
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z > 4 (e.g. ISteidel etTil Il998t IPentericci et"ail l200d: 
IVenemans et all l2002t IShimasaku et all l2003|) "The dis- 
covery of these proto-clusters with sizes of several Mpc 
was made possible by deep wide-field optical surveys, 
concentrating on a narrow redshift interval usi n g the 
Lyman-break or narr ow-band Lva llKurk et all I2OO0I: 
IVenemans et all I2OO2I) or Ha l|Kurk et all l2004|) excess 
techniques, followed by multi-slit spectroscopy on 8- 10m 
telescopes. Sim ilar surveys have also been done in 
dom fields (e.(^ . lHu et allll99a[l2nn2HRhoads et all 



ran- 



JR 

n. l200l . 



iKodaira et ail l2003HShimasaku et al.l l200 

Hierarchical galax y formation models {e.g. 
iKauffmann et all Il999f) predict that the best fields 
to search for such overdensities are those containing a 
massive galaxy. High redshift radio galaxies (HzRGs; 
z > 2) are therefore ideal targets, as they are not only 
among the most massive galaxies known at high z (e.g . 
iDe Breuck et all l2002t iRocca-Volmeranee et all \20o4 - 
but their lower redshift (z > 0.5) cou nterparts are 
also l ocate d in c l uster environment s {e.g. iHill fc Lillvl 
Il99lt iBestL l200nt iBest et all kOtt^ . This motivated an 
intensive search of Lya emitting companions around six 
radio galaxies with 2.1 < z < 4.1, leading to the discovery 
of overdensities of Lvct emitters of 5 to 15 compared with 
the random fields l|Venemans et all I2OO2I Ebo3il . 

However, the Lya excess technique detects only a 
fraction of the companion objects. For example, in the 
proto-cluster of Lyma n break galax i es at z=3.09 in the 
SSA 22 region (LBG; ISteidel et all Il99?j) . only 25% of 
the LBGs have a sufficiently high Lya eq uivalent width 
to be included in a Lya excess selection l|Steidel et all 
This incompleteness is also illustrated by the pres- 
ence of a similar overdensity of Ha emitters and ex- 
tremely red objects in the proto-cluster surrounding the 
z=2.16 radio galaxy PKS 1138—262, while these galax- 
ies are mostly not detec ted using the Lya excess tech- 
nique l|Kurk et all l20o3) . Deep Chandra observations of 
this field have also found an overdensity of X-ray sources, 
of which at least 2 ar e AGN with i n the p roto-cluster 
jPentericci et all l2002() . ISmail et all l)2003b(i also report 
the detection of four Chandra X-ray sources coincident 
with submm sources surrounding three HzRGs. 

These proto-clusters also contain a large amount of 
gas, as revealed by the Lya haloes surrounding the 
HzRGs, which ha ve physical scales up to > 200 kpc (for a 
recent review, see lvan Breugel et alll200^ . Together with 
the increased merg e r rate s in higher redshift clusters {e.g. 
Ivan Dokkum et all Il999l) . this provides the ingredients 
to induce wide-spread starbursts and AGN, which could 
be revealed by their thermal dust (sub-)mm emission. 
Statistical overdensities of (sub-)mm galaxies (SMGs) 
have indeed been found from SCUBA bolometer imag- 
ing of the fields surrou nding the z=3.09 'redshift spike' 
f Chapman ct al.','200l'), the z=2.39 radio galaxy 53W002 
(■Smail et al. 2003a), the z=3.8 ra dio galaxy 4C 41.17 



2=1.8 QSO llStevens et all I2OO4I) 
HzRG fields {siivs^mSrESoJ"^ 



I Ivison et all 
and six other 
this paper, we present a 1.2 mm map covering the 



central 25 arcmin^ of the most distant proto-cluster 
known to date, surrounding the z=4.1 radio galaxy 
TN J1338-1942. We find an overdensity of 1.2 mm 
sources, which we identify with optically very faint galax- 
ies using a deep 1.4 GHz map, but find no overlap between 
the population of excess Lya emitters and the 1.2 mm 
sources. Throughout this paper, we use a A— cosmology 



with 

3rge1 et aLl 
minosity distance is -Dl = 
to 7.0 kpc. 



nM=0 .27 
II2Q0I 



and f7A=0.73 
At z=4.1,the lu- 



=37.65 Gpc, and 1" corresponds 



2. Observations and data reduction 

2.1. I RAM 1.2mm imaging 

To image the field of TN J1338-1942 at mm wavelengths, 
we used the 37- and 117-channel M ax Planck Bolometer 
arrays (MAMBO-1 and MAMBO-2; iKrevsaet all Il998|) 
at the IRAM 30m telescope on Pico Veleta, Spain. 
MAMBO has a half-power spectral bandwidth from 210 
to 290 GHz, with an effective bandwidth centre for steep 
thermal spectra of ^250 GHz (1.2 mm). The effective 
beam FWHM is 10'.' 7 with an array size of 4'. 

The observations were done in a pooled observing 
mode during the winter 2001-2002 season. Due to the low 
declination, the field could only be observed for 4 hours in 
a given night, with elevations between 29° and 33°. The at- 
mospheric zenith opacities at 1.2 mm varied between 0.12 
and 0.25. The total on source integration time was 17.0 
hours, of which 2.6 h were obtained with the 117-channel 
array, and 14.6 h with the 37-channel array. We used the 
standard on-the-ffy mapping technique, comprised of 41 
subscans of 40 s each, while chopping the secondary mir- 
ror in azimuth at 2 Hz. To minimize the residual effects of 
the double beam point spread function, we used different 
chop (wobbler) throws (39", 42" or 45") and/or scan direc- 
tions for each map. We used 7 different pointings offset by 
80" to cover the entire field of the VLT/F0RS2 imaging to 
uniform depth. We checked the pointing and focus at least 
once per hour using the bright point source 1334—127, and 
found the pointing to be stable to within <2". The abso- 
lute flux calibration is based on observations of several 
standard calibration sources, including planets, resulting 
in an estimated accuracy of 15%. 

\ ye ana l yzed t he data using the MOPSI software pack- 
age (|Zvlkalll998l) . We subtracted the skynoise, and com- 
bined the double-beam maps using a shift-and-add proce- 
dure, producing for each map a positive image bracketed 
by two negative images of half the intensity located one 
chop throw away. Because we combined our 37 maps ob- 
tained with different chop throws, the effect of the con- 
fusion due to negative sidelobes is minimized, but still 
present in regions of high source density. The noise level 
increases outward in our co-added map of the field. In 
Fig- n two contours show the region within which the rms 
noise level is less than 1.2 and 0.6 mJy (before smoothing), 
enclosing areas of 25.6 and 2.6 square arcmin, respectively. 
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RfGHT ASCENSION (J200Q} 

Fig. 1. VLT i?— band image fgrevscale: IVenemans et 311120021) with MAMBO 1.2 mm signal-to-noise map (smoothed 
to 11") overlaid as thin/blue contours. Contour levels are -3, -2, 2, 3, 4 and 5a, with a the local rms noise level 
(negative contours are dashed). The two large purple contours delineate the regions with a <0.6 and 1.2 mJy/beam. 
The MAMBO beam size is indicated in the lower left corner. The red open crosses indicate spectroscopically confirmed 
Lya companions to the radio galaxy TN J1338— 1942, and the green open cross labeled Rl an additional FR II radio 



2.2. VLA radio imaging 

To obtain accurate positions of the mm sources, and to 
search for possible radio- loud AGN counterparts, we ob- 
served the T N J1338-1942 field with the Ver y Large 
Array (VLA; iNapier. Thompson fc Eke^ Il983^ on UT 
2002 April 1 to 12 for a total of 12 hours in the A-array 
at 20 cm. We observed in a pseudo-continuum, spectral 
line mode with 7 x 3.125 MHz channels. Wc monitored the 
point source 1351 — 148 every 40 min to provide amplitude, 
phase and bandpass calibration, and used an observation 
of 3C 286 to provide the absolute flux calibration. 



We performed standard spectral-line calibration and 
editing of the data using the NRAO AIT'S package, 
and employed standar d wide field imaging techniques 
l)Tavlor. Carilh. fc Perle v. 1999). The final 7'.5x7'.5 image 
has an rms noise level of 15 /iJy beam^^, except in the 
area close to the central radio galaxy, which is limited by 
the ability to clean the bright radio source (see Fig. 
The FWHM resolution of the restoring beam is 2'.'3 x 1'.'3 
at a position angle PA=0°. 
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Table 1. Astrometry of the MAMBO, VLA and VLT sources in the TN J1338-1942 field. 



MAMBO VLA VLT 

Source RA(J2000) DEC(J2000) RA(J2000) DEC(J2000) RA(J2000) DEC(J2000) 

h m s o I It h m s o I II h m s o I II 



MOl 


13 


38 


17.56 


— 19 


42 


30.7 


13 


38 


17.48 


— 19 


42 


29 


5 


13 38 


17.43 


— 19 


42 


29.9 


M02a 


13 


38 


25.65 


—19 


41 


25.5 


13 


38 


25.47 


— 19 


41 


21 


5 


13 38 


25.45 


— 19 


41 


22.1 


M02b 














13 


38 


25.37 


— 19 


41 


26 


5 












M02c 














13 


38 


25.81 


-19 


41 


26 


2 












M03 


13 


38 


25.97 


-19 


42 


33.9 


13 


38 


26.10 


-19 


42 


31 


3 


13 38 


26.06 


-19 


42 


30.7 


M04 


13 


38 


26.84 


-19 


42 


22.8 


13 


38 


26.83 


-19 


42 


25 


8 


13 38 


26.82 


-19 


42 


26.2 


M05 


13 


38 


37.79 


-19 


41 


48.8 


13 


38 


37.91 


-19 


41 


49 


1 


13 38 


37.92 


-19 


41 


48.9 


M06 


13 


38 


30.13 


-19 


41 


06.5 


13 


38 


30.04 


-19 


41 


04 


5 


13 38 


30.06 


-19 


41 


04.7 


M07 


13 


38 


29.89 


-19 


41 


36.9 


























M08 


13 


38 


28.86 


-19 


43 


20.9 
















13 38 


28.77 


-19 


43 


27.5 


M09 


13 


38 


28.33 


-19 


42 


13.2 


























MIO 


13 


38 


29.78 


-19 


39 


32.1 


13 


38 


29.81 


-19 


39 


31 


9 


a 




a 






Rl 














13 


38 


22.77 


-19 


42 


30 





13 38 


22.76 


-19 


42 


29.9 



" This source is located outside the VLT image boundary. 



Table 2. Photometry of the MAMBO and VLA sources in the TN J1338-1942 field. 



Source 


S'20cm 


<S'l200Mm 




5'450fim 


i?[3"] 


NB[S"] 


K[3"] 






mjy/beam 


mjy 


mJy 


mag 


mag 


mag 


MOl 


148±34 


6.2±1.2 


10.1±1.3 


21.5±6.4 


26.8±0.3 


>28 


a 


M02a 


54±25 


4.1±0.8 


6.6±1.3'' 


12.5±6.0'' 


25.6±0.1 


26.7±0.7 


a 


M02b 


46±25 








>28 


>28 


a 


M02c 


47±26 








>28 


>28 


a 


M03 


120800±4300 


2.3±0.6 


6.1±1.3'^ 


25.3±9.3 


22.38±0.01 


19.2±0.1 


19.97±0.05 


M04 


47±17 


2.0±0.5 






24.08±0.03 


24.9±0.2 


18.72±0.03 


M05 


91±25 


3.8±0.8 


9.9±1.2'' 


14.9±5.9* 


25.06±0.07 


25.6±0.3 


a 


M06 


41±24 


2.3±0.6 


5.7±1.3 


<20 (3cr) 


27.8±0.8 


>28 


a 


M07 


<50 (3(t) 


4.0±0.6 


3.3±2.0 


<33 (3a) 


>28 


>28 


a 


M08 


<50 (3cr) 


2.4±0.7 






25.26±0.09 


26.0±0.4 


19.2±0.1 


M09 


<50 (3cr) 


2.3±0.6 






d 


d 


a 


MIO 


195±28 


3.2±1.0 






a 


a 


a 


Rl 


10900±400 


<1.8 (3cr) 






25.6±0.1 


25.6±0.3 


a 


Source located outside the 


coverage of the VLT/Keck images. 









SCUBA flux densities could be underestimated due to the extent of the 1.2mm source (see §3.2 and Fig. I^L 

We quote the S/N weighted average of our value (Sgsoum = 3.1 ± 1.8) and the value quoted bv iReuland et all (|20^ 
(5850A.m =6.9±1.1). 
No clear optical/near— IR counterpart. 



2.3. JCMT submm photometry 

We obtained 850 fim and 450 fim photometry of 6 
MAMBO sources previously identified using the VLA and 
VLT imaging (§3.1) using the Submillimetre Comm on- 
User Bolometer Array (SCUBA; iHoUand et all. Il999(l on 
the 15m James Clerk Maxwell Telescope on UT 2003 
February 17 to 22, for a total integration time of 16 hours 
(including overheads) spread over 6 sources. The respec- 
tive beamsizes are 14'.'7 at 850 /im and 7'.'5 at 450 //m. We 
used the recommended 9-point jiggle photometry mode, 
which samples a 3 x 3 grid with 2" spacing between grid 
points, while chopping in azimuth by 60". The atmo- 
spheric conditions were excellent, with zenith opacities at 
850 fim varying between 0.11 and 0.24. We checked the 
pointing and focus roughly once per hour using the bright 
point source 1334—127, and found the pointing to be sta- 



ble within <3". The absolute flux calibration is based on 
observations of CRL618 and Mars. 

We reduced the data in the stan dard way usiii g 
SCUBA User Reduction Facility (SURF: I.Tennesl Il997l) . 

We removed individual data samples that were deviant 
from the mean by >Aa from the data set, flatfielded, cor- 
rected for extinction, despiked, and subtracted the resid- 
ual sky level from the off-source bolometers. Finally, we 
calculated the mean and rms from the time sequence to 
remove all points deviating by more than 3f7 from the 
mean. We adopt the mean of the remaining data as our 
final estimate of the intensity. 

2.4. Optical and near— I R imaging and spectroscopy 

Optical and near-IR imaging data of the TN J1338 — 1942 
field exist in the literature. iDe Breuck et al.l l|l999i) pub- 
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lished a NIRC/Keck iiT— band image c overing the cen- 
tral I 'xl', centred on the radio galaxy. IVenemans et alJ 
l|2002l) published deep i?— band and narrow-band Lya 
images obtained with the FOcal Reducer/low dispersion 
Spectrograph (FORS) at the ESO Very Large Telescope 
(VLT), covering 6.'4x6f2, and also spectroscopically con- 
firmed 20 Lya emitters in this field, of which 15 fall 
within the cr < 1.2 mJy region of our MAMBO map. 
Oyerzier et al. (in prep.) also obtained a ii'— band im- 
age of the south-eastern corner of the proto-cluster with 
the Infrared Spectrometer And Array Camera (ISAAC) at 
the VLT. In the following, we shall use both the iiT— band 
and i?— band images to identify the optical counterparts of 
the mm and radio sources. We corrected all magnitudes 
for Gal a.ctic extinction E(B — T/)=0.097 us ing the dust 
maps of lSchlegeL Finkbeiner. fc Davisl lll998l) and th e ex- 
tinction cunfeofJCarfeli^^l^^ l|l989l) . 

Optical spectroscopy of three sources (MOl, M05 and 
M08) was obtained as part of two VLT multi-slit spec- 
troscopy programs. Details about these observations will 
be given in future papers (De Breuck et al., in preparation; 
Overzier et al., in preparation). 

3. Results 

3.1. Extraction of the mm sources 

We use the MAMBO signal-to-noise map (Fig. to 
identify candidate mm sources within the region where 
a <1.2 mJy/beam. To assess the reliability of the candi- 
date sources, we also made 'probability maps' by splitting 
up the data sets, and looking for gross discrepancies be- 
tween them. This would identify spurious sources appear- 
ing in only part of the data. The reliability of our positive 
detections is further strengthened by the absence of neg- 
ative peaks < — 4tT in the S/N map (dashed contours in 
Fig.^. The three most significant negative sources are all 
located in between M02 and M07, consistent with them 
being residual double beaming effects (see §2.1). 

Table ^ lists the positions of 10 candidate 1.2 mm 
sources found in the MAMBO map, and Table [3 their 
flux densities. We fitted the sources in the map smoothed 
to 11" to Gaussians unconstrained in width, and in cases 
where the fitted width is smaller than the beam size, to 
constrained width Gaussians of 11" FWHM. The peak 
fluxes quoted in Table[21are averages of Gauss fits in maps 
smoothed to 11" and 12", and to eliminate a base problem 
due to double beams surrounding the brighter sources, in 
a map truncated toward negative values at 0. The quoted 
uncertainties are quadratic sums of the dispersion in these 
three fits and the rms noise level at the position of the 
source in the 11" resolution map. 

3.2. Identification of the mm sources 

Due to the large MAMBO beam (FWHM ^ 11"), we 
often find several potential optical counterparts within 
the positional uncertainty. Higher resolution interferomet- 



ric imagin g would be needed for more unique identifica - 
tions {e.g. iDannerbauer et alll2002tlDunlop et alll20o'^ . 
Comparisons of positions derived from MAMBO maps and 
from the VLA or IRAM Plateau de Bure interferometer 
yiel d median offsets of 2-3" with a tail extendi ng to 7- 
8" l|Eales et all 120031: IDannerbauer et all 12004 - To ob- 
tain more accurate positions, we can take advantage of 
the tight rad io-to- far- infrar ed correlation in star-forming 
galaxies {e.g. Cgn^Qg , 1^2^ ), which appears to hold out to 
high redshift l|GarrettLl20fl2^ . Because the faint radio syn- 



chrotron emission traces the same region as the mm emis- 
sion, we can use deep sub-arcsecond resolution 1.4 GHz 
VLA maps to piiipoint the dust emission to within <1" 
{e.g. Ilvison et allEioi. 

The surface density o f radio sour c es br ighter than 
30 ^Jy is ~2.5 arcmin"^ ijlvison et al.L l20o3) . so we ex- 
pect to find only 0.035 such sources within a radius of 4" 
from the MAMBO positions. Because 30fi3y is at the 2a 
noise level in our radio map, we expect to find 0.4 positive 
2a peaks within 4" from the nominal MAMBO positions. 
In an attempt to discriminate against noise peaks, we also 
checked for optical or ncar-IR counterparts for the ra- 
dio sources. We consider an optical/near-IR source with 
a >2a radio peak within the astrometric uncertainty of 
0"4 as a likely counterpart to a MAMBO source, if it lies 
within 4" from the MAMBO position. In the following, we 
discuss the individual identifications of the 10 candidate 
MAMBO sources, and of an additional radio source in the 
field. Table ^ lists the mm, radio and optical positions, 
and Table 13 the photometric measurements. 
MOl (Fig.EJ: This is the brightest source in our MAMBO 
map. The identification is secure with a 4.4cr radio source 
coinciding with an i?=26.8 source, which we used as the 
basis for SCUBA photometry. We have attempted deep 
VLT/FORSl spectroscopy on this object (De Breuck et al. 
in preparation), but did not detect any emission, so the 
spectrum is not dominated by strong line emission from 
an AGN. 

M02 (Fig. 0): This MAMBO source appears to consist 
of two components; an unconstrained Gaussian fit yields 
a size of 19"7xl7"0. Within a radius of 5" from the 
MAMBO position, there are three 1.4 GHz sources with 
flux densities >5'i.4ghz >45 /iJy, i.e. at signiflcance levels 
of ~2f7 (see Table 121 , while we expect to find only one 
such source. This suggests a complex, merging structure 
of galaxies, as hinted by the diffuse i?— band source coin- 
cident with the brightest of the three VLA sources, M02a 
(see Tables n and 121 . We have used this position for the 
SCUBA photometry, but the relatively low Ss5o / 51200 ra- 
tio indicates that the SCUBA beam may have missed part 
of the fiux, suggesting M02a is not the correct identifica- 
tion of the brightest 1.2 mm emission. Only interferomet- 
ric (sub-)mm imaging could provide a better understand- 
ing of this interesting source. 

M03 (Fig.|2): This 3.8cr MAMBO sourc e lies 3"2 from the 
host g alaxy of the powerful radio source ijDe Breuck et all 



'^1999'), and is therefore most likely related. We used the 
if— band position to obtain MAMBO on-off observations 
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Fig. 2. Keck/NIRC /T— band image (greyscales) with the MAMBO 1.2 mm signal-to-noise map (smoothed to 11") 
overlaid as thick/blue contours and the VLA 1.4 GHz map as thin/red contours. Contour levels for the MAMBO map 
are 2, 3 and 4a, with a the local rms noise level. Contour levels for the VLA map are -0.1275, -0.09, -0.063, -0.045, 
0.045, 0.063, 0.72, 4.08, 5.775, 62.28 and 130.56 mJy/beam. The MAMBO and VLA beam sizes are indicated in the 
lower left and lower right corners, respectively. The open cross indicates the position of the optical/radio identification. 



and SCUBA photometry, but both of these gave lower 
than expected flux densities (see Table EJ, suggesting that 
the dust emission could be offset from the iiT— band posi- 
tion, as has been seen in the z=3.79 radio galaxy 4C 60.07 
llPaoadoDOulos et all I2OOOI) . In Table El we quote the 
S/N weighted average of our 5's.';n//m=3.1zl il.8 mJy mea- 
surement and the S'85o^tm=6.9±l.l mJy of iReuland et al.1 
ll20Q4l). We do not averag e the 450 fim flux densities, as 
the iReuland et all l|2004|) photometry was obtained un- 
der adverse atmospheric conditions leading to a nominal 
measurement of — —36.3 ± 31.5 mJy. Such dis- 

crepant Ss50fj.m flux densities between observations at dif- 
ferent epochs have been noticed before for MG 2141-1-192 



llArchibald et all l200lt IReuland et all l2004l) _aiid_ sug- 
gest submm variability. But, as IReuland et al.l l)2004l) ar- 
gue, the extended nature of the submm emission excludes 
this explanation, unless there is a significant contribu- 
tion from non-thermal AGN emission. However, a power- 
law extrapolat ion from the total 8.2 GHz flux density 
l)Pentericci et al.. 2000) predicts a negligible AGN syn- 
chrotron contribution of only 30 /iJy at 250 GHz, but given 
the unusually asy r nmetr ic morphology of the radio source 
toe Breuck et al.lll999|) . we cannot exclude the presence 
of a variable, flat spectrum radio core contributing to the 
submm flux densities, as has been seen in B2 0902+343 
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Fig. 3. VLT i?— band image (greyscales) with the MAMBO 1.2 mm signal-to-noise map (smoothed to 11") overlaid as 
thick/blue contours and the VLA 1.4 GHz map as thin/red contours. Contour levels for the MAMBO map are 2, 3, 4, 
5, 6 and 7a, with a the local rms noise level. Contour levels for the VLA map are -42,-30, 30, 42, 60 and 85 /iJy/beam. 
The VLA beam size is indicated in the lower right corner. The open cross indicates the position of the optical/radio 
identification. 



ijPownes et al.llT996|) . High resolution imaging at frequen- 
cies of several tens of GHz would be needed to test this. 
M04 (Fig. 121): This 4cr MAMBO source lies S'.'O from a 
faint radio source coincident with a very red galaxy with 
R — iir=5.36±0.04, which is a likely identification. 
M05 (Fig. 01: We obtained SCUBA photometry of the 
extended radio and optical source, confirming the identi- 
fication. This is the only source which has spatially ex- 
tended 1.4 GHz emission; the close correspondence be- 
tween the diffuse optical and radio morphologies strongly 
suggests the 1.4 GHz emission traces starburst rather 
than AGN emission. Deep VLT/FORSl spectroscopy (De 
Breuck et al. in preparation) also revealed a faint feature- 
less continuum, but no emission lines. An unconstrained 



Gaussian fit to the MAMBO S/N map yields a size of 
20'.'8xl8'.'0, but the S/N of our detection is only 3.7 in a 
11" smoothed map, so the detection of spatially extended 
emission is tentative at best. 

M06 (Fig. El : We obtained SCUBA photometry at the 
position of a very faint radio and optical identification, 
confirming the reality of this faint MAMBO source. Note 
that the apparent spatial extent of the MAMBO emission 
cannot be trusted because this source is only detected at 
the 3.4(7 level. 

MOT (Fig. 13): Although this is one of the brightest 
MAMBO sources, there is no obvious radio or opti- 
cal identification. We have obtained SCUBA photom- 
etry of the radio/optical source at RA=13''38"30.22^ 
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Fig. 3. , continued 

DEC=-19°41'33'.'06, but obtained only a 1.5a signal of 
5'850/im=3.7±2.5 mJy, while at the nominal MAMBO po- 
sition, we obtain a slightly higher signal (see Table[2Jl. We 
do not have a good candidate optical or radio counterpart 
for this MAMBO source. 

M08 (Fig. OJ: This faint MAMBO source lies 6'.'7 from 
an extremely red object (ERO) with iir=19.2±0.1 and 
R - is:=6.1±0.15. We have obtained deep F0RS2/MXU 
spectroscopy of this ERO (Overzier et al., in preparation), 
detecting a faint red continuum and an emission line at 
8I20A, which we tentatively identify as [Oil] A 3727 at 
z=1.18. Note that from blank-field ERO surveys, the den- 
sitY of objects with i?—j<'> 6 and K<19.2 is 0.10 arcmin"^ 
i Daddi et all l200(il) . so the chance of finding such an ob- 



ject within 6'.'7 is P=OA%. It is therefore possible that 


M08 is a dusty ERO ( e.q. ICimatti et al 




1998: Dcv et all 


Il999l ISmail et all l2002bl: Takata et al. 




2003). However. 



because the ERO is not the closest possible identification, 
it is statistically not the most likely identification. Note 
that the apparent spatial extent of the MAMBO emission 
cannot be trusted because this source is only detected at 
the 4(7 level. This may also indicate that this source has 
made it into our sample due to the confusion of two sources 
too close to be detected individually. Deeper mm/submm 
observations would be needed to verify this, and to deter- 
mine if the dust emission is related to the ERO or not. 

M09 (Fig. 121: We find no radio source within one 
MAMBO beamsize, and no clear optical/near— IR iden- 
tifications. 

MIO (Fig. EI): This MAMBO source falls just outside the 
VLT _R— band image. It coincides with a strong 0.2 mJy 
radio source. 

Rl (Fig. fH) : This is a moderately bright radio source with 
a FR II iFanaroff & RilevL \l974) morphology We detect 
no 1.2 mm emission from this source. 
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Fig. 4. MIO: MAMBO 1.2 mm signal-to-noise map 
(thick/blue contours; smoothed to 11") and the VLA 
1.4 GHz map (thin/red contours). Contour levels for the 
MAMBO map are 2 and 3a, with a the local rms noise 
level. Contour levels for the VLA map are -42,-30, 30, 42, 
60, 85, 120 and 170 ^Jy/beam. The VLA beam size is in- 
dicated in the lower right corner. This source falls outside 
the VLT _R— band image. 

3.3. Extended 1.2mm emission 

At least one (M02) of the three sources detected with 
S/N>5, and possibly the weaker source M05 appear to 
have significantly extended 1.2 mm emission. M02 appears 
to consist of two barely resolved components, though none 
of them have a clear candidate optical or radio identifi- 
cation. Spatially extended submm emission has been re- 
ported before in several HzRGs and at least two compan- 
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ion sources lllvison et al l l2000HSteveiis et allEool Ivison 
et al., in prep.). This suggests that star formation in these 
objects occurs over scales of several tens of kpc. 

A possible alternative explanation is that M02 (and 
maybe M05) could consist of multiple ima ges by stron; 
gravi t ational lensin g by a foreground cluster l|Kneib et al 
I2nn4t iBorvs et a1.L \2004\ . However, we have not de 
tected any gravita tional arcs in the ACS/ HST images 
l)Milev et al.r l2004j) or in the VLT image (Fig. P). We 
therefore consider it extremely unlikely that gravitational 
lensing can explain the spatial extent in the mm emission. 
A more detailed analysis of the lensing properties, includ- 
ing an analysis of a possible weak shear using the ACS 
image will be presented in a future paper (Overzier et al., 
in prep.). 

3.4. Surface density 

Although many of the MAMBO sources were detected at 
low S/N level, most of them appear to be real because 
because they have been confirmed with pointed sub-mm 
photometry, or because they have plausible counterparts 
at radio and/or optical wavelengths. We detect 10 sources 
brighter than 2.0 mJy (peak S/N>3, a level at which we 
should be close to being complete). Omitting M03, which 
corresponds to the pre-selected radio galaxy in the field, 
this corresponds to a surface density of 0.35 arcmin~^ 
(roughly 1% of the confusion limit). Note that this value 
is likely to be an overestimate because we consider sources 
down to the 3a level , whic h may have lead to flux boost- 
ing {e.g. lEales et al.U2003|) . This effect raises intrinsically 
fainter sources above the detection threshold due to the 
addition of noise (instrumental and atmospherical) or due 
to confusion of sources too faint to be detected individu- 
ally. The la t ter aff ect may be happening in M08 (Fig.|31). 
lEales et all l|2003|) argue that for MAMBO surveys, the 
average boosting factor is 14%, which is slightly lower 
than for the SCUBA surveys, possibly due to the smaller 
MAMBO beam. 

We now compare our surface density with the inte- 
grated 1. 2 mm source counts fr om the MAMBO blank field 
surveys l|Bertoldi et al.U2000|) . Using the 1.2 mm counts 
from Bertoldi (private communication) and ICreve et alJ 
lj2004|) . our 5*1.211101 source density (without flux boost- 
ing correction) is roughly twice as high as expected from 
the random-field density. However, if we consider only 
the 3 sources with 5i.2nini >4 mJy, the source den- 
sity is 0.12 arcmin^^, three times higher than expected. 
Assuming a Poissonian distribution, this represents a 7a 
overdensity. These bright sources are all detected with a 
local^ S/N>5, so they will not be affected by flux boost- 
ing. This overdensity suggests that most of the brighter 
MAMBO sources and a few of the fainter ones may be 
related to the proto-cluster. Such statistical overdensities 



Table 3. Photometric redshift estimates. 



Source 


2est from 


2;ost from S'850Mm/'S'l200Mm 


MOl 


^•^-0.9 




M02 


2.8tll 


7 1 +00 
' -^-4.9 


MOS" 




7 1+" 

' — 1.9 


M04 


2.2+2} 


M05 


n 9+1.8 
^■^-1.0 


9 1 +6.8 
^•^-2.1 


M06 


9 9+4.1 
^■^-1.0 


2 4+" 
^•^-2.4 


M07 


>2.9 


>4.2 


M08'' 


>2.3 




M09 


>2.3 




MIO 


-■-•^-0.7 





^ Note that the uncertainties in Table |5| also include the fit- 
ting uncertainties (see §3.1), and are hence larger than the local 
rms. 



Radio galaxy at 2:=4.1. 

Possibly identified with a 2spcc=1.18 ERO, see §3.2. 



of SMGs a re also seen in SCUB A maps of other HzRG 
fields {e.g. IStevens et all l2003h . In the next section, we 
use our radio, mm and submm photometry to constrain 
the redshifts of these MAMBO sources. 



3.5. Photometric redshift estimates 

We attempted very deep (tint=14h) optical spectroscopy 
of MOl and M05, but did not detect features to determine 
the redshifts. None of the MAMBO sources show excess 
narrow-band emission in the 4.075 < z < 4.123 coverage 
of the Lya filter, and only M04 and M08 have /iT-band 
photometry, so we cannot use conventional optical photo- 
metric redshift techniques. As this is a common situation 
for submm galaxies, several redshift estimators have been 
de veloped based on the submm and radio data only. 

ICarilh fc YunI l)l999|) proposed to use the radio-to- 
submm spectral index as a redshift estimator, which 
is based on the local radio-to-far-infrared correlation. 
Table El hsts the red shift estimates calculated using the 
ICarilli fc Yu"^ l|2000|) aj^^^^^ models. The quoted uncer- 
tainties include both the la measurement uncertainties 
in our photometry and the la scatter in the models. 
The estimated redshifts are clearly lower than the z=4.1 
of the radio galaxy, though z—A.l is still within the la 
uncertainties for most sources (except MIO). However, 
a comparison with spectroscopic redshifts of SMGs in- 
dicates that the spectral-index redshifts tend to be sys- 
tematically underesti mated for ^an pr > 2 objects (e.g . 
IClements et"all l2004j) . In fact, the Icarilli fc YunI l|200CI) 
models predict a spectral index a^^4Qy| = 1.00 for z=4.1, 
so for the average 5*1200 =3.3 mJy in our maps, we expect 
to find a radio source with a fiux density 5i.4GHz = 18A<Jy, 
well below the dete ction limit of our VLA map. However, 
IPetric et all l)2003l) report 1.4 GHz detections of ~70^Jy 
in two z > 5 radio quiet quasars with 5i2oo flux den- 
sities of 0.9 and 5.5 mJy. A possible explanation for this 
higher than expected 1.4 GHz emission is that an optically 
undetected AGN contributes to the radio emission. Note 
that AGNs have been spectroscopi cally confirmed i n two 
proto-clustcrs surrounding Hz RGs l |Le^evre et al.l Il996t 
iPentericci et al.. .2002,'! , while ISmail et al.i (,2003b|l report 
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the discovery of Chandra X-ray sources coincident with 
submm sources surrounding three HzRGs, suggesting the 
MAMBO sources surrounding TN J1338-1942 ~ if they 
are part of the protocluster - may weU contain AGNs. 
Hence, a relatively bright radio detection of a MAMBO 
source does not exclude it as a potential member of the 
z = A.l protocluster. 

lEales et all ()2003|) have also predicted the redshift 
evolution of the S'850Mm/'S'i200Atm ratio, by fitting a two- 
temperature model to a sample of 104 galaxies from the 
IRAS bright galaxy survey iDunnc & Ealcs, 2001). We 
have used their median predicted value (Fig. 4 in the 
lEales et al.l l2003l paper) to obtain an additional redshift 
estimate. Table |21 lists these estimates; the quoted uncer- 
tainties include the la measurement uncertainties in the 
5'850Aim/'S'i200/jm ratios and the full spread in the model 
predictions. We find that, except for a z > 2.2 limit for 
MOl and M02, this ratio does not provide a useful con- 
straint on the redshifts due to (i) the small difference in 
wavelength between the mm and submm points and (ii) 
the relatively low S/N of our detections. 

In summary, while the radio, mm and submm photom- 
etry are generally consistent with z=4.1, the uncertainties 
from these redshift estimate techniques are far too large 
to provide proof that the sources are at the redshift of 
the radio galaxy. Better sampled SEDs, especially in the 
submm would be needed to constrain the redshifts with 
photometry only. 

3.6. Statistical 1.2 mm and 1.4 GHz emission from the 
Lya emitters and Mean Star Formation Rates 

None of the 14 spectroscopically confirmed Lya emit- 
ters within the a < 1.2 mjy region (excluding the ra- 
dio galaxy) is detected at > 2a in our 1.2 mm map. To 
test if there is a statistical signal in our non-detections, 
we have stacked the emission from these 14 positions, 
and find (S'i.2mm)=0.25±0.24 mJy. Using this 3a up- 
per limit, and assuming dust parame t ers T ri=33 K and 
/3=2.0 used for LBGs jBaker et all l200lj) . we derive 
(^fir) ^ 4.4 X IO^^Lq, implying a mean star fo rmation 
rate (SFRfir) < 440 Mgyr"! {e.g. lOmont et all EoOl). 

Similarly, from our VLA map, we find 
(S'i.4GHz)=— 0.2±2.9 /iJy. Using the relation between 
the radio emission and star formation rate ( Condorj, 
Il992l) . we can use this 3a upper limit of 8.7 /iJy on 
the 1.4 GHz emission to calculate an upper limit to 
the star formation rate (SFR) from the Lya emitters. 
Assuming a radio spectral index aradio=— 0.8, we derive 
(SFRradio) <265 Moyr-i. 

The values should be compared with the SFR derived 
from the mean Lya luminosity of the same 14 emitters 
(excluding the radio galaxy) (LLya)=2.77x lO^^erg s~^. 
Assuming a case B ratio Lya/Ha=1 0, and using the 
iKennicutt. Tamblvn. fc CongdonI l)l994j) relation between 
SFR and Lh^, we find (SFRLya)=2.2 Mgyr"!. Note that 
this value is likely to be an underestimate, as Lya is often 



quenched by dust emission, as illustrated bv iKurk et alJ 
(2004), who report Lya/Ha ratios significantly lower than 
the case B value for a sample of Lya and Ha emit- 
ters surrounding the z—2.16 radio galaxy PKS 1138—262. 
However, it is obvious that the deep Lya imaging probes 
much lower SFR than the MAMBO and VLA maps. 

4. Discussion 

Although we could not put strong constraints on 
the redshifts of the 9 MAMBO sources surrounding 
TN J1338— 1942 using photometric redshift estimators or 
very deep VLT/FORSl spectroscopy of two of them, the 
analysis of the source density and the photometric redshift 
estimates suggests half of the 9 MAMBO sources, and in 
particular the brightest ones, may well belong to the z=4.1 
proto-cluster. To confirm (or refute) this requires alter- 
native redshift determinations such as (i) deep near-IR 
spectroscopy, (ii) mid-IR spectroscopy with IRS /Spitzer 
using the PAH features, and (iii) mm spectroscopy using 
molecular CO lines. 

Figure ^ shows that the nine MAMBO sources are 
not distributed uniformly around the radio galaxy, al- 
though the sensitivity of our MAMBO map radially de- 
creases from the radio galaxy. The four brightest MAMBO 
sources (MOl, M02, M07 and M05) are all north of the 
radio galaxy, while the densest area of Lya emitters is 
locat e d southeast of the radio galaxy ijVenemans et all 
I2OO2I: iMilev et all l2004|l . This suggests that the cen- 
tre of the proto-cluster is not necessarily to the south- 
east of the radio galaxy, as suggested by the distribu- 
tion of Lya emitters. None of the 14 spectroscopically 
confirmed Lya emitters in the MAMBO field were de- 
tected at 1.2 mm or 1.4 GHz. Similarly, none of the 
MAMBO sources show excess Lya emission (if they are 
within the 4.075 < z < 4.123 coverage of the narrow-band 
Lya filter). A possible explanation would be that their 
large amounts of dust, as traced by the 1.2 mm emission 
may have quenched the Lya emission (e.g. iKunth et all 
^98), putting them below the narrow-band excess cutoff 
(EWicst >15 A). This apparent absence of overlap be- 
tween both populations shows the importance of using 
multiple wavelength techniqu es to obtain a mo r e com plete 
picture of t h e pro to-cluster. IChapman et alJ l|2000() and 
IWebb et alJ (|2P0^ also report a <20% overlap between 
SMGs and LBGs, wi th a few n otable excep t ions (e.g. 
IChaoman etali |2Qo3)- However, IWebb et alJ l(2003l) do 
find a large cross-clustering correlation amplitude between 
both populations in the largest of their 3 fields. We do not 
detect enough MAMBO sources in the TN J1338-1942 
field to perform such a cross-correlation analysis, but our 
detection of a statistical overdensity of MAMBO sources 
does suggest some r elation to the confirin ed overden- 
sity of Lyg emitters (" Venemans et al.l l2002|) and LBGs 
llMilev et all [2004 ) in the same field. 

The faintness of the Lya in our objects also con- 
trasts the publis hed spectra of SMGs, which often show 
bright Lya lines IjChapman et aL. .2003,) . However, these 
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SMGs are on average at lower redshifts (z) ^ 2.4 
ijChapman et all l2003f) . while our MAMBO sources are 
potential members of the z=4.1 proto-cluster. Although 
we only have K—hand information for 3 sources, which are 
all relatively bright class -II (I — K < 5, K <21) or class- 
I (/ -- K > 5) sources iflvison et al.l l200otlsmail et all 
l2002a|) . 3 of the 4 brightest MAMBO sources may well 
be very faint class-0 {K > 21) sour ces which are much 
harder to obtain redshifts from {e.g. IPannerbauer et aD . 
I20n2t iFraver et all l2004|) . either because they are more 
highly obscured, or because they are at higher redshifts. 
Our success in finding higher redshift objects may have 
been helped by the use of a longer selection wavelength 
(1.2 mm instead of 850 /zm). Indeed, based on the low ra- 
tio of 850 fim to 1.2 mm flux, several authors argue that a 
significant fraction of mm galaxies may be at z > 3 {e.g. 
lEales et all 120031: lAretxaga et all l2003l). Most MAMBO 
sources in the NTT deep field |Dannerbauer et aD . l2002(l 
are also class-0 sources, supporting the trend that 1.2 mm 
selected sources appear to be fainter (and hence maybe at 
higher redshift) than 850 /J,m selected sources. 

The detection of X-ray emission in several bright 
submm sources wit hin the proto-cluste rs surrounding 
three other HzRGs (jSmail et al.l l2003b|) suggests that 
AGN may be present within those sources. Half of the 
spectra presented by IChaoman et al.l l|2003^ also show 
type-II AGN lines, but we failed to detect these in our 
very deep VLT spectroscopy of MOl and M05. Hence, 
we have no indication that AGN are present in our 
MAMBO sources, but their emission l i nes rn av well have 
been obscured {e.g. iReuland et al.L l2003() . Assuming 
the MAMBO sources are at z=4.1, and dust parame- 
ters Td =50 K and (3=1.5 l)Benford et alill999|) . we de- 
rive LpiR ~ 2 X IO^'^Lq X 5'i.2mm/mJy in the range 
4-12xlO^^L0, implying star formation rates of 4000 to 
12000 Mf7^vr ~^. if the dust is entirely heated bv s t ar for - 
mation {e.g. lOmont et allElOl: De Breuck et alll2003l) . 
For M08, which is possibly at z=1.18, we derive Lfir ~ 
1.6 X lO^^L© and a star formation rate of 1600 MQyr~^. 
Such hi gh star formation rates have been reporte d for 
HzRGs ((Archibald et al.[l200ltlR,euland et alll2004l) . but 
our MAMBO map indicates that they may also occur out 
to distances as far as 2 Mpc from the central radio galaxy. 
If the MAMBO sources are really at z=4.1, this would 
suggest that these proto - cluste rs have multiple, possibly 
aligned (jPentericci et al.[ l2002fl concentrations. However, 
the VLT narrow-band Lya image shows that none of the 
mm sources seen in the TN J1338— 1942 proto-cluster have 
huge Lyg haloes, like those seen around the c e ntral HzRG 
{e.g. IVenemans et all l2002t iReuland et all l2003(l . The 
HzRG therefore appears the best candidates to evolve into 
the present-day giant elliptical. 

5. Conclusions 

We summarize the results from our multi-wavelength 
observations of the z=4.1 proto-cluster surrounding 
TN J1338-1942 as follows: 



• We detect 10 candidate mm sources with peak fluxes 
having S/N>3 in our MAMBO map. Of these, at least 
eight sources with 5*1.211101 >1.3 mJy have possible ra- 
dio and/or optical/near-IR counterparts, and 5 are con- 
firmed at S/N>4 with pointed SCUBA submm photome- 
try. Three sources have 51.2111111 >4.0 mJy, while comparing 
with source counts from blank field surveys, we expect to 
find only 1 such source in the unassociated field popula- 
tion. 

• The radio-to-submm and mm-submm photometric 
redshift estimates do not provide strong constraints on 
the possible redshifts of the MAMBO sources, although 
for 9 of the 10 sources, they are consistent with z=4.1 
within the uncertainties. 

• One of the faint MAMBO sources is possibly related 
to an ERO with R — K=6.1, which has a likely spectro- 
scopic redshift of z=1.18, implying a star formation rate 
up to 1600 Moyr-^ 

• None of the 14 spectroscopically confirmed Lya emit- 
ters show detectable millimetre emission, and supposing 
they are at z=4.1, none of the 1.2 mm sources show an ex- 
cess Lya emission, indicating no apparent overlap between 
both populations. 

• The mean star formation rate of the 14 spectroscop- 
ically confirmed Lya emitters, as derived from the deep 
Lya imaging is two orders of magnitudes lower than the 
upper limits derived from the stacked VLA and MAMBO 
maps, illustrating that the radio and mm maps probe 
much more actively star- forming galaxies. 

Our multi- wavelength observations of this z=4.1 
proto-cluster suggest that the Lya excess technique does 
not detect the massive starburst companions within 
2 Mpc. Their high SFR could be triggered by interactions 
between different proto-cluster members, as indicated by 
the diffuse nature of the optical identifications of M05 and 
M02. 
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